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We present an sp 3 tight-binding model for the calculation of the electronic and optical properties 
of wurtzite semiconductor quantum dots (QDs). The tight- binding model takes into account strain, 
piezoelectricity, spin-orbit coupling and crystal-field splitting. Excitonic absorption spectra are 
calculated using the configuration interaction scheme. We study the electronic and optical properties 
of InN/GaN QDs and their dependence on structural properties, crystal-field splitting, and spin- 
orbit coupling. 

PACS numbers: 78.67.Hc, 73.22.Dj, 71.35.-y 
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I. INTRODUCTION 



Group-Ill nitrides have shown great potential in op- 
toelectronic devices with a wide range of applications. 1 
InGaN/GaN quantum wells typically constitute the ac- 
tive region in light-emitting diodes and laser diodes, cov- 
ering a wide spectral range from near ultraviolet to in- 
frared. In particular, devices with an active region based 
on pure or almost pure InN are of great interest to reach 
operational frequencies in the infrared spectral ranged 
The small InN band gap (0.7-0.8 eV)£ can be extremely 
useful for telecommunication- wavelength devices. The 
combination of this inherent property of InN with the 
self-assembly of nanostructures based on this material, 
provides further possibilities for useful future optoelec- 
tronic devices. In particular, quantum dot (QD) struc- 
tures are promising candidates, acting as electron-hole 
recombination centers increasing the emission efficiency. 
These zero-dimensional nanostructures also have the po- 
tential to act as single-photon emitting devices 4 

This work is dedicated to the investigation of the elec- 
tronic and optical properties of self- assembled InN QDs. 
The present study is based on a fully atomistic em- 
pirical tight-binding model. In contrast to multi-band 
k • p approaches it takes into account the structure of 
the underlying atomic lattice and is capable to describe 
the electronic wave functions beyond an envelope func- 
tion approximation. Our model includes strain effects 
and electrostatic built-in fields, which are of major im- 
portance in group-Ill nitride based nanostructures with 
an underlying wurtzite crystal structure. In this pa- 
per, special attention is paid to the possible influence 
of the weak crystal-field (CF) splitting and weak spin- 
orbit (SO) coupling on the localized single-particle wave 
functions as well as on the optical properties of the inves- 
tigated structures. These effects have commonly been ne- 
glected in previous studies, some of which were dedicated 
to GaN/AIN nanostructures^^, others to InN/GaN 
nanostrucutrcs 9, 10, x - i 12 i 13 i 14 . Since in these materials 



both effects are merely of the order of a few meV (Refj^), 
this approximation can be well justified. 

On a more fundamental level the resulting symmetry 
properties dictated by spatial and spin degrees of freedom 
and, taking SO coupling into account, by the coupled in- 
fluence of both, have the potential to support or lift cer- 
tain degeneracies in the electronic energy spectra. Based 
on group-theoretical arguments we discuss that in the 
system under investigation for both electrons and holes, 
at most two-fold degeneracies are supported. This is es- 
pecially of interest with regard to a recent discussion of 
the energy level structure in semiconductor QDs with a 
zinc blende structure^ and recent results for wurtzite In- 
GaN/GaN QDs obtained within an 8-band kp model 16 . 
Following these general group-theoretical arguments, our 
numerical results show that SO coupling and CF splitting 
slightly change the results for the InN/GaN QD system 
under investigation. However, the additional splittings 
in the electronic energy shell structure caused by the SO 
coupling are at most of the order of a few meV. From this, 
no significant qualitative changes of the optical proper- 
ties are found: The excitonic absorption lines show only 
very small additional splittings and the interband dipolc 
selection rules are basically unaltered. 



II. THEORY AND QUANTUM DOT MODEL 

A. The Tight-Binding Model 

For the investigation of the single particle states in 
small QD structures, the description by a multi-band 
approach is required. In order to take into account the 
underlying wurtzite structure of the structures under in- 
vestigation, we choose a microscopic sp 3 tight-binding 
(TB) model. The general aspects of the TB model are 
discussed in detail in Ref™,. Here, we briefly summarize 
the main ingredients of this model, and focus our atten- 
tion on the role of crystal-field splitting and spin-orbit 
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FIG. 1: Schematic band structure of wurtzite semiconductors 
with conduction band and three valence bands. The valence 
band splittings introduced by crystal field splitting and spin- 
orbit coupling are denoted by Ai and A2, respectively. Ai 
and A2 are calculated according to Eq. ((TJ. Additionally, the 
symmetries of the irreducible representations Fi are given. 



coupling. 

In the sp 3 TB model of RefJ^ the relevant electronic 
structure of anions and cations is, for each spin orien- 
tation, described by the outermost valence orbitals, s, 
p x , p y and p z , and the overlap of these basis orbitals is 
restricted to nearest neighbors. Being only of the order 
of a few meV, so far, the influence of CF splitting and 
SO coupling has been neglected in the model. In the 
present work, we extend our TB models to introduce 
these two contributions and investigate their possible in- 
fluence on the electronic single-particle states and ener- 
gies in InN/GaN QD systems. 

The combination of CF splitting and SO interaction 
leads to a so-called three-edge structure in the vicinity 
of the r point. The top of the resulting valence band 
structure is commonly labeled as A, B, and C bands in 
order of increasing energy. This three-edge structure is 
schematically shown in Fig. [TJ Two of these three bands 
are of and one of Tg symmetry. To describe this band 
structure in the vicinity of the T point in the framework 
of a TB model one has to take into account both CF 
splitting and SO interaction. 

The SO coupling is included as outlined by Chadi in 
Re£i&. Due to the high ionicity of the bonds in the ni- 
tride system-^, the contribution of the SO coupling to the 
valence band structure is dominated by the anion contri- 
butions. Therefore, by introducing the parameter A, we 
include spin-orbit coupling at the anion sites only. The 
parameter A is used to reproduce the correct splitting Ai 
of the valence bands A and B. 

As discussed in Ref.—, the small CF splitting A c f of 
the wurtzite crystal differentiates the p z orbital from the 
p x and p v orbitals. Pseudo-potential calculations in lo- 
cal density approximation indicate that for the studied 



TABLE I: Tight-binding parameters (in eV) for the nearest 
neighbors of wurtzite InN and GaN. The notation of Ref.— is 
used. 
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materials the bulk crystal field splitting between the A 
and C valence bands, schematically shown in Fig. [2j 
cannot be reproduced from first principles, unless third- 
nearest- neighbor interactions are taken into account 20 . 
The TB model discussed in Ref considers only nearest- 
neighbor hopping matrix elements and treats the four 
nearest neighbor atoms as equivalent. To account for the 
A — C splitting within the empirical sp 3 TB model with 
nearest-neighbor coupling, we introduce the additional 
parameter p ^ on the anion sites for the on-site ma- 
trix elements of the p z orbitals. This additional term is 
used to reproduce the splitting A2 of the A — C bands at 
the zone center T. 

With four atoms per unit cell, the resulting Hamilto- 
nian is a 32 x 32 matrix for each k-point. This Hamil- 
tonian parametrically depends on the different TB ma- 
trix elements. The parameters given in Ref. 12 have been 
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FIG. 2: (Color online) Bulk band structure of wurtzite (a) 
InN and (b) GaN obtained using an sp 3 TB model including 
crystal field splitting and spin-orbit coupling. The insets show 
the three-edge valence band structure in the vicinity of the V 
point. 
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re-calculated to reproduce the three-edge structure the 
wurtzite band structur o 3 ' 21 ! 22 in the vicinity of the T 
point. The additional matrix elements A and pz are 
adjusted to reproduce the splittings between the different 
valence bands (A, B and C bands), which are given by 23 

(i) 

where SO and CF splitting are denoted as A so and A c f 
respectively. Table U summarizes the resulting TB pa- 
rameters. The bulk band structures obtained from these 
parameters are shown in Fig. [2] for GaN and InN. The 
'complicated' valence band structure in the vicinity of 
the r point shows very good agreement with other TB 
models 24 (cf. insets of Fig. (2), ab-initio approaches 25 , 
pseudo potential and k • p calculations^. 



B. Geometry of the Quantum Dot Structure 

In the framework of a TB model, the QD is modeled 
on an atomistic level. The TB parameters at each atom 
site R of the underlying wurtzite lattice are set according 
to the bulk values of the respective occupying atom. In 
the following, we consider lens-shaped InN QDs, grown 
in c direction and residing on an InN wetting layer (WL). 
The entire structure is embedded inside a GaN matrix. 
For the numerical calculations, a finite cell (box) with 
a wurtzite lattice and fixed (zero) boundary conditions 
is used. The cell is sufficiently large to avoid numerical 
artifacts in the bound single particle states due to the 
cell boundaries (in particular, no artifacts from the arti- 
ficial cubic symmetry of the supercell are found) . In the 
following we discuss three different QD sizes with diame- 
ters d = 4.5, 5.7, 7.7 nm and heights h = 1.6, 2.3, 3.0 nm, 
respectively. The WL thickness is one lattice constant c. 

In a first step we neglect strain induced displacements 
of the atoms and concentrate on effects which can ex- 
clusively be attributed to crystal-field splitting and spin- 
orbit coupling. The influence of strain effects will be 
discussed separately in Sec. IIIIBI 

In contrast to semiconductor heterostructures with a 
zinc blende structure, the III-V wurtzite nitrides exhibit 
a considerably larger built-in electrostatic field 2 ^. In or- 
der to account for this field in our model, the electrostatic 
potential 4> p (r) is obtained from the solution of the Pois- 
son equation and enters as a site-diagonal contribution 
V p (r) = —e<fi(r) to the TB Hamiltonian 1 ^. The polar- 
ization P has two contributions, the spontaneous polar- 
ization p s P ont and the piezo-electric contribution pP lczo 
caused by strain inside the system. For the latter we 
apply the approximation described in Ref. 2 - and assume 
ppiezo ^ e ^ ^ w hich is a reasonable approximation for 
the considered QD geometry^. Further details of this 
procedure are given in Refji 2 -. 



C. Many-Body Hamiltonian, Coulomb and Dipole 
Matrix Elements 

Having discussed the TB Hamiltonian used for the cal- 
culation of the bound single-particle states, we now turn 
our attention to the investigation of the optical properties 
of the studied QD system. We start with the following 
Hamiltonian H that describes the dynamics of the inter- 
acting charge carriers in the system: 

H = H a + H c + H D . (2) 

This Hamiltonian consists of three parts and is given in 
the basis of the QD one-particle eigenstates. The contri- 
bution Hq 

i i 

is the one-particle part, which is diagonal in the chosen 
basis, 

H c = ^V%u44c k c l + ~^V£Hh\h]h k ht 

ijkl ijkl 

-^V^uhtcXckfH, (3) 

ijkl 

describes the Coulomb interaction of electrons (e) in the 
conduction band states and holes (h) in the valence band 
states, and 

H D = '£(e(i\Er\j)c i h j +h.c.) , (4) 

includes the coupling of the electronic system to an exter- 
nal electromagnetic field E in dipole approximation. The 
creation and annihilation operators for electrons (holes) 
in the single-particle state \i) with energy e\ (ej 1 ) are de- 
noted by c\ (h\) and c i (hj), respectively. The Coulomb 
interaction matrix elements are labeled by V$jy • 

The calculation of the Coulomb interaction matrix el- 
ements requires - at least in principle - the knowledge 
of the localized basis states implicitly underlying the TB 
wave functions. However, since the Coulomb matrix ele- 
ments are dominated by the long-range character of the 
interaction, in the calculation of these matrix elements 
the charge densities in the localized orbitals are approx- 
imated by point charges. A more detailed discussion of 
this issue is given in Ref* 1 ^. This approximation leads 
to the following explicit form of the Coulomb matrix el- 
ements: 

Vw = Err- E a p c^^^ViK - R) , (5) 
with 

F(R-R0 = 47reoEr f R -R^| fOT R^R' 

and 

V(0) = =4- /" d 3 rd 3 r'- d e \ rj « V . 
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The expansion coefficients c l a R are related to the i th one- 
particle wave function $,(r) = ^ Q R c l a R </> a ,R.(r) where 
Qj R,(r) denotes the atomic wave functions localized at 
the lattice site R. 



To calculate the dipole matrix elements d 



ell 



3C 



(ipf entering Eq. (j4|), we use numerically orthog- 

onalized Slater orbitals 1 ^ to account for the short-range 
character of the dipole operator. The numerically or- 
thogonalized Slater orbitals fulfill the basic properties of 
the localized basis states underlying the TB model: sym- 
metry, spatial orientation 2 ^, and orthogonality. We also 
include the anion-cation structure of the crystal and the 
slight nonlocality of the dipole operator by including con- 
tributions from up to second nearest neighbors. 

In analogy to the bulk systems, a separation of the 
orbital and spin part (both included in the index a) is 
prohibited by the spin-orbit coupling. In contrast to the 
bulk will be discussed in the following section, the 

strong band mixing prevents a strict classification of QD 
single-particle states according to their angular momen- 
tum. Therefore, even total angular momentum selection 
rules are no longer applicable. Ignoring the band mix- 
ing characteristics in the zero-dimensional structures, any 
treatment of many-body effects based on strict selection 
rules for the total angular momentum, yields inaccurate 
predictions of level degeneracies. However, as discussed 
in detail in RefJ^, the selection rules can always be ana- 
lyzed on symmetry grounds. 

With the dipole and Coulomb matrix elements and the 
many-body Hamiltonian, given in Eq. @, the calcula- 
tion of optical spectra can be carried out as described in 
RefsM. 



III. RESULTS FOR LENS-SHAPED InN 
QUANTUM DOTS 



Elections 
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A. Single Particle states and Energies 

Having determined the TB parameters, the single- 
particle states and energies of the three different QDs 
discussed in Sec. Ill Bl can be calculated. 

First we turn our attention to the single particle states 
of the large QD (d — 7.7 nm, h — 3.0 nm). In order to 
assess the impact of the CF splitting and SO coupling, 
we have performed our calculations in three steps. In a 
first step we neglect both SO interaction and CF effects 
(A so = 0, A c f = 0). In this case we are left with the sp 3 
TB model discussed in RefJ^. In step two we introduce 
only the crystal field splitting (A so = 0, A c f ^ 0), by in- 
cluding the additional parameter „ in the TB model. 
In the final step both contributions are taken into account 
(A so ^ 0, A c f ^ 0). In each step the TB parameters are 
re-optimized in such a way that band gap and the ener- 
getic positions of other bands of the wurtzite bulk band 
structure at the T point are reproduced. In other words, 
in the first step (A so = A c f = 0) we use the parameters 



FIG. 3: (Color online) Top view of the large lens-shaped InN 
QD structure with the first three bound states for electrons 
(upper part) and holes (lower part). Depicted are isosurfaces 
of the probability density with 10% (blue) and 60% (red) of 
the maximum value. 



given in Ref. 12 . In step two (A so = 0,A c f ^ 0) the pa- 
rameter Ep z p is included in the TB-model to reproduce 
the splitting of the A — C bands. Of course all the other 
matrix elements are also re-adjusted to obtain the correct 
band structure. Starting from these TB parameters, in 
the final step (A so ^ 0, A c f ^ 0) the additional parame- 
ter A is taken into account to reproduce the splitting of 
the bands A — B. Again, all parameters are re-adjusted 
to reproduce the correct energetic positions of the differ- 
ent bands. The resulting parameters are given in Tab. [T] 
Figure [3] shows the QD geometry and first three bound 
one-particle states for electrons and holes, respectively, 
including the influence of the built-in field. 

According to the nodal structure, the depicted electron 
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ground state ipf can be classified as s-like. The first two 
excited states and can be classified as p + and p_ 
states, respectively. Such a classification is not possible 
for the hole states, since these states undergo strong band 
mixing effects. Considering only a single valence band for 
the description of the bound hole states in an InN QD is 
not valid. The observation of strong band mixing effects 
is in agreement with other multi-band approaches ! 5 ' 16 ! 31 
From Fig. [3] we can deduce that CF splitting and SO 
coupling do not alter the single-particle level structure. 
In other words, for the ordering of the first three bound 
electron and hole states in a lens-shaped InN QD, the 
contributions from SO and CF splitting are negligible. 

After this discussion of the single-particle states we fo- 
cus on the single-particle energies of the large QD. In 
Tab. IIIII the energies of the first three bound electron 
and hole states under the influence of the built-in field 
are displayed. From this table we conclude that the elec- 
tron states are only slightly shifted to lower energies by 
CF splitting and SO coupling. Without SO coupling 
(A so = 0), and taking only CF splitting into account 
(A c f ^ 0), the hole states are shifted to higher ener- 
gies. With SO coupling, the hole energy spectrum is 
shifted to lower energies, compared to the case neglect- 
ing both contributions (A c f = 0, A so = 0). Additionally 
it turns out that the degeneracy of the hole states 
and ips (p-shell) is lifted when SO coupling is considered. 
Of course, each state is still twofold degenerate due to 
time reversal symmetry 3 - 3 ". Because of the small SO en- 
ergies of the bulk materials, the splitting is rather small 
(A^h — 1.2 meV). The same is true for the electron 
p-states. Here the influence of the SO coupling is even 
weaker (A^h^h =0.1 meV). 

As discussed recently^, for the system under investi- 
gation, a TB model which neglects the weak crystal-field 
splitting and spin-orbit coupling, must result in degener- 
ate p-shells for electrons and holes. The origin of these 
degeneracies is the C% v symmetry of the combined system 
of QD geometry and underlying wurtzite lattice. 

From the splitting of the hole states ^2 and 4>3 , one can 
deduce that the spin-orbit interaction alters the symme- 
try of the system. This can be understood by an analysis 
of the corresponding character tables for C^ v single and 
double groups. 

Without spin-orbit coupling, the symmetry of the sys- 
tem is determined by the single group C% v . Looking at 
the character table, shown in Tab. [IIJ we conclude that 
this group allows for double degenerate levels, since the 
group contains a two-dimensional representation T3. One 



TABLE II: Character table for the single group C 3 „ (Ref.— ) . 



{E} {2Cs} {3a v } 





1 


1 


1 


r 2 


1 


1 


-1 


r. 


2 


-1 






TABLE III: Single-particle energies for the large InN QD in 
the presence and absence of crystal field splitting and spin- 
orbit coupling. Each of the given states is two-fold degenerate 
due to time reversal symmetry. 





A so = 0, A cf = 


A so = 0, A c( ^ 


A so ± 0, A cf ^ 


Ef [eV] 


1.4770 


1.4585 


1.4557 


Ei [eV] 


1.6660 


1.6464 


1.6417 


El [eV] 


1.6660 


1.6464 


1.6418 


E\ [eV] 


0.9021 


0.9041 


0.8993 


£2 [eV] 


0.9021 


0.9041 


0.8981 


E\ [eV] 


0.8964 


0.8989 


0.8924 



example of such degenerate states are the p-shell states 
for electrons (V>f; V 7 !) and holes (ipi ; ^2 X shown in Fig. [3] 
Of course, in the absence of the SO coupling, in addition, 
each state is two-fold spin degenerate. 

Including SO coupling, one has to deal with the dou- 
ble group Cs v . The character table of the double group 
is given in Tab. IIVI This group allows only two dimen- 
sional representations, even if the time reversal symmetry 
is included^ 3 -. More specifically, the degeneracy of the ir- 
reducible representations L 3 and L 4 is not doubled by 
the time reversal symmetry— In other words, no four- 
fold degenerate state in the energy spectrum can exist. 
Consequently, the electron p states (^f and V'f) are also 
no longer degenerate, but in this particular case the split- 
ting is only of the order of some fj,eV. We note that this 
is in contrast to the findings in Ref, 16 , where (at least 
numerically) exact degeneracy of the electronic p-states 
has been reported. 

A central result of the previous wor k 12 ' 13 was that the 
strong internal electrostatic field can reverse the ener- 
getic ordering of the first three bound hole states. We 
find that, for the intermediate and the largest InN QD, 
in the presence of the built-in field, the ground state is 
formed by the twofold degenerate p states ip^ and , 
shown in Fig. [3] This behavior is interchanged with de- 
creasing QD size, where, for the smallest QD, the s state 
■03 becomes the hole ground state. To concentrate on 
this reordering of the hole s and p shell, the energy split- 
ting ABjj, = E' s l - Ep is displayed in Fig. [3J In order to 
analyze the impact of CF splitting and SO coupling on 
the ordering of the hole level structure, we calculate the 



TABLE IV: Character table for the double group C 3v (Ref.- 24 ). 



{E} {E} {2C 3 } {217 3 } {3a u } {3a,} 



r a 


1 


1 


1 


1 


1 


1 


r 2 


1 


1 


1 


1 


-1 


-1 


r 3 


2 


2 


-1 


-1 








r 4 


2 


-2 


1 


-1 








r 5 


1 


-1 


-1 


1 


i 


— i 


L, 


1 


-1 


-1 


1 


—i 


i 



6 



• -No CFSO 
■-CFSO 

♦ - CFSO+Strain 




FIG. 4: 
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(Color online) The energy splitting AE* 
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the absence (No CFSO) and presence (CFSO) of crystal field 
(A c f) and spin-orbit (A so ) splitting. Results are also shown, 
where in addition to crystal field and spin-orbit splitting a 
strain field is included (CFSO+Strain). The influence of the 
strain field on the electronic properties is discussed in detail 
in Sec. IIIIBI In all three cases, AE$ p changes sign with 
increasing QD diameter d, as a level reordering occurs (the 
dashed lines are included as a guide to the eye) . The built-in 
field is included in all the calculations. 



splitting AE^p when both contributions are introduced 
in the TB approach. Since the hole p shell is no longer 
degenerate, we average over the single-particle energies of 
the states Vt and tp!j: . The energy splitting AE^ p with 
SO coupling and CF splitting is also shown in Fig. 2] 
From the comparison with AE^ p in the absence of these 
contributions, we find that the SO coupling and the CF 
splitting have only a negligible effect on the energy split- 
ting AE^ p . Furthermore, the ordering of the first three 
bound hole states is unaffected by SO coupling and CF 
splitting. 

In summary, the CF splitting alone cannot alter the 
symmetry of the system and leads only to a small energy 
shift of the first three bound electron and hole states. The 
SO interaction, and only this contribution, can modify 
the symmetry, and lifts certain degeneracies. However, 
the splitting of the electron and hole p shell due to the SO 
coupling, is very small compared to the level spacing of 
the different shells. Moreover, in the presence of SO cou- 
pling and CF splitting, one obtains the same level order- 
ing of the energetically lowest hole states (s- and p-shell) 
as in the case were these contributions are not taken into 
account. This analysis indicates that it is well justified 
to neglect these small corrections of the CF splitting and 
the SO coupling in the system under consideration. 



B. Influence of strain 

So far we have neglected the influence of the lattice 
mismatch between InN and GaN. The lattice mismatch 



TABLE V: Material parameters for GaN and InN. If not in- 
dicated otherwise, all parameters are taken from Ref.—. 



Parameter 


GaN 


InN 


a (A) 


3.189 


3.545 


c(A) 


5.185 


5.703 


A so (eV) 


0.010 (Ref£) 


0.005 (Ref£) 


A cf (eV) 


0.017 (Ref4) 


0.040 (Ref£) 


Cm (GPa) 


11.4 


9.4 


C 33 (GPa) 


38.1 


20.0 


(a c -.Di) (eV) 


-9.6 (Ref. 36 ) 


-9.6 (Ref.-S 6 -) 


(oe - Da) (eV) 


-8.2 (Ref. 36 ) 


-8.2 (ReiM) 


D 3 (eV) 


1.9 (Ref.^) 


1.9 (Ref^) 


D 4 (eV) 


-1.0 (Ref .22) 


-1.0 (RefM) 


a c (eV) 


-4.9 (Ref£) 


-3.5 (Ref£) 



leads to the appearance of a strain field in the nano- 
structure. This field modifies the energies of the bound 
electron and hole states. 

For the electron states, caused by the underlying 
wurtzite lattice and the assumed QD geometry, the strain 
field produces only an energy shift of the bound single 
particle states^. In contrast to QDs with a zinc blende 
structure^, no degeneracies are lifted by the strain field. 
The situation is more complicated for the hole states. 
The strain field modifies the local valence band edges, 
and can therefore lead to a splitting of the different en- 
ergy bands. As discussed in RefJ£, due to the biaxial 
strain in the basal plane, the first two valence bands (A 
and B) are shifted to higher energies, compared to the 
unstrained material, whereas the third valence band (C 
band) is shifted to lower energies. These energy shifts 
may also increase the influence of spin-orbit coupling and 
crystal field splitting on the bound single particle states. 

To investigate the influence of these shifts and the ef- 
fect of the possible valence band splittings, we proceed 
in the following way: 

Since the TB parameters are fitted to the bulk band 
structure, we re-calculate the bulk band structure of InN 
in the presence of a strain field. To obtain the strain 
dependent valence band edge we apply 37 



E- 



2.3 



A cf - 

A cf 
2 



Aso 

3 



e f + K 



(6) 
(7) 



\ 



+ ^(A S o) 2 , (8) 



where e and X e are given by 

6e = D 3 e zz +D 4 {e x 
K = D 1 e zz +D 2 {e x 



+ € yy) 
+ e vv) 



Here, D±, D2, D3 and D4 are the valence band deforma- 
tion potentials and en denotes the diagonal components 
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FIG. 5: (Color online) Calculated TB valence band structure 
in the vicinity of the F point. Unstrained (black solid lines) 
and under compressive biaxial strain in the basal plane with 
elastic relaxation along the [0001] direction (red dashed lines). 
The splitting between the A and B valence bands at the V 
point is w 3.1 meV (« 3.1 meV) in unstrained (strained) InN. 



of the strain tensor e. For the components m we assume: 

2C 13 



ao — a 



a 



C: 



:S:S 



(9) 



where ao and a are the lattice constants of the substrate 
and the QD material, respectively, and C13 and C33 are 
the stiffness constants. Strictly speaking, these equations 
apply only to a quantum well, since they neglect shear 
strain components and the fact that the strain field in a 
QD is position dependent. However, due to the symme- 
try of the QD structure and the underlying wurtzitc lat- 
tice, the strain field can not lift degeneracies. This is also 
confirmed by results of Winkelnkemper et. al^ for an 
In^Gai-^N QD with a comparable symmetry. Further- 
more, the hole states, for which the modification of the 
valence band edge is of major importance, are strongly 
localized in the region of the QD. Therefore, the variation 
of the strain field outside the QD is of minor importance. 
In the region of the InN QD, we use the lattice constant 
a of GaN, to take into account, that the QD is pseudo- 
morphically grown on the GaN substrate. 

According to Refill, the band gap shift is given by 



AE sap 



Pr, 



with 



Pc, 



= a c e z 



- a c (e x 
D 2 (e. 



(0« + A e ) 



+ £yy) , 

: + e yy) ' 



(10) 



where a c denotes the conduction band deformation po- 
tential. The different parameters are listed in Tab. [V] 
The TB parameters are readjusted to reproduce the bulk 
band structure in the vicinity of the T point. The result- 
ing valence band structure in the vicinity of the T point 
is depicted in Fig. [5l The band structure exactly reflects 
the behavior which is expected for the local band struc- 
ture in the region of the QD: 16 the first two valence bands 



TABLE VI: Single-particle energies for the large InN QD in 
the presence and absence of strain effects. Each of the given 
states is two-fold degenerate due to time reversal symmetry. 
The internal electrostatic field is included in the calculation. 





A so / 0, A cf / 




Without Strain With Strain 


Ef [eV] 


1.4557 


1.7740 


El [eV] 


1.6417 


1.9625 


El [eV] 


1.6418 


1.9626 


E\ [eV] 


0.8993 


0.9126 


El [eV] 


0.8981 


0.9118 


E^ [eV] 


0.8924 


0.9069 



(A and B) are shifted to higher energies while the third 
band (C) is shifted to lower energies. The calculated TB 
parameters are now used to investigate the influence of 
strain effects 011 the electronic states. 

Starting from these new TB parameters one can re- 
calculate the single-particle states and energies of the 
three different QDs discussed in the preceding section, 
taking into account both CF and SO splitting. First 
we focus on the single-particle energies of the large QD 
(d = 7.7 nm, h = 3.0 nm). In Tab. lVIl the energies of the 
first three bound electron and hole states, including the 
built-in field are displayed. Each of the given states is 
two-fold degenerate due to time reversal symmetry. The 
results in the absence and in the presence of the strain 
effects are compared. First of all, the strain merely pro- 
duces an energy shift of the single particle states. The 
strain field shifts both electron states and hole states to 
higher energies. This behavior reflects the local band 
edge shifts of conduction and valence bands. Following 
the discussion of the previous section, no four- fold degen- 
erate states can exist taking SO coupling into account. 
As already discussed, the first two excited states i?| and 
i?| are nearly degenerate in the absence of the strain field. 
Also in the presence of the strain field, these states are 
split by the SO coupling by 0.1 meV only. This analysis 
shows that the splitting of the electron p-states is not al- 
tered by the strain field and remains very small compared 
to the energy separation of electron s- {ipf) and p-shell 
(1P2 and f^l). The splitting of the hole states i/'i and ^2 
is nearly unaffected by the strain field. 

Additionally, the level ordering is not modified by the 
strain field. In the presence of the strain field, the intrin- 
sic electrostatic field still reverses the energetic ordering 
of the first three bound hole states. The energy splitting 
A^ p between the hole s and p shell including the strain 
field is displayed in Fig. 0J We compare these results 
with the other results shown in Fig. 2J where the strain 
field is absent. In conclusion, the strain field is of minor 
importance for the energy splittings, and, most impor- 
tantly, the ordering of the first three bound hole states is 
unaffected. 
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FIG. 6: (Color online) Excitonic absorption spectra for the 
largest InN QD in the absence of crystal field and spin-orbit 
splitting (a), in the presence of crystal field splitting and ab- 
sence of spin-orbit coupling (b) and in the presence of crystal- 
field and spin-orbit splitting (c). The insets show the same 
data for frequencies close to the absorption peaks in (b) and 
(c). 



and iplj: is very small, the splitting of the two peaks on 
the low energy side is also very small. This emphasizes 
again, that in the system under consideration the SO 
coupling and the CF splitting introduce only negligible 
corrections to the excitonic spectrum. Due to the sym- 
metry of the QD and the underlying wurtzite structure, 
there is no polarization anisotropy. This is in contrast 
to lens-shaped InAs QDs with a zinc blende structure, 
as discussed in RefJ^, where the Ci v symmetry leads to 
polarization anisotropy. 

The peak on the high energy side mainly corresponds 
to the excitation of the hole in the state tp^ and the 
states ip2 an d "03 ■ Since the electron states i/j-z an d ^3 
are exactly degenerate in cases (a) and (b) and nearly 
degenerate in case (c), only a single line is visible on the 
high energy side (cf. inset of Fig. [6J. 

Including strain effects, we obtain nearly the same 
splitting of the states V'i and ■ Furthermore, the elec- 
tron states 1P2 an( i V'f are still nearly degenerate. Only 
the single-particle energy gap is enlarged by the strain 
field and therefore the whole excitonic absorption spec- 
trum is shifted to higher energies. The excitonic absorp- 
tion spectrum (not shown) resembles the spectrum in 
Fig. [6] (c). 



IV. CONCLUSION 



C. Excitonic absorption spectra 

In Sec. Ill Cl we have discussed the calculation of dipole 
and Coulomb matrix elements. The evaluation of exci- 
tonic absorption spectra in this section can be performed 
starting from the many-particle Hamiltonian, Eq. ([2|), in 
second quantization as given in Sec. Ill CI For the lo- 
calized states configuration-interaction calculations are 
performed. For the sake of simplicity, only the first three 
bound electron and hole states are included. This can 
be justified by their energy separation to higher shells 
in the structure. The excitonic absorption spectra are 
calculated using Fermi's golden rule.— 

Figure [5] (a) shows the excitonic absorption spectrum 
for the largest InN QD in the absence of CF and SO split- 
ting, while Fig. [6] (b) displays the spectrum in presence 
of the CF splitting but in the absence of SO coupling. 
In Fig. [5] (c), the absorption spectrum in the presence 
of both CF and SO splitting is depicted. The different 
absorption lines in each spectrum correspond to the ex- 
citation of an exciton in the QD. In case (a) and (b) the 
peak on the low energy side corresponds to transitions 
where the electron is mainly in the ground state tpf, and 
the hole is mainly in the states ■f/'i an d ip2- Since the 
hole states and are degenerate one obtains only a 
single peak on the low energy side. Due to the SO cou- 
pling, the states ^1 and i/^ are split by about 1.2 meV. 
This splitting in the single particle states results in two 
lines on the low energy side in Fig. [6] (c) (see also the 
inset). Since the splitting in the single-particle states ip^ 



In this work we have presented an atomistic TB cal- 
culation of the electronic and optical properties of lens- 
shaped InN QDs. We focused our attention on the influ- 
ence of the crystal-field splitting and the spin-orbit cou- 
pling on the electronic structure as well as on the optical 
properties. As it turns out, only the spin-orbit coupling 
lifts certain degeneracies in the single particle spectrum. 
This result is confirmed by the inspection of the charac- 
ter table of the double group C^v , which reveals that no 
four fold degenerate state can exist. However, from our 
calculations we obtain only small splittings in the single- 
particle spectrum as well as in the excitonic absorption 
spectrum. Our results indicate that it is a reasonable as- 
sumption to neglect spin-orbit coupling and crystal field 
splitting in case of lens-shaped InN QDs with a wurtzite 
structure. 

We do not rule out that the importance of spin- 
orbit coupling and crystal-field splitting may depend on 
the specific system under investigation even within the 
group-Ill nitride material system where both effects are 
intrinsically weak. 
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